            figure; axesm('globe')

%            gridm('GLineStyle','-','Gcolor',[.8 .7 .6],'Galtitude', .02)

           load coast
            plot3m(lat,long,.01,'k','LineWidth',2)
            view(3)
              base = zeros(180,360); baseref = [1 90 0];


      hs = meshm(base,baseref,size(base));
      colormap gray
      camlight right
      material([.8 .9 .4])
   % material([.5 .5 .5])
set(gca,'Box','off')
hidem(gca)
view(60,30)
[a,b] = xlsread('c:\manoj\projects\ace\stations.xls');
plotm(lat,long,.01,'k','LineWidth',2)
plotm(a(:,1),a(:,2),'r.','markerSize',30)
load c:\manoj\geomag\dip_eq.mat;
plotm(lat111,long111,'b','LineWidth',2);